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Stress  Failure  Distribution 


The  failure  of  an  element  of  equipment  such  as  a  structural  member 
or  a  soldered  joint  due  to  overloading  is  in  general  of  the  nature  of 
chance  failure  (Ref.  1),  since  (neglecting  long-term  effects  such  as 
fatigue)  if  the  beam  survives  the  stressing,  it  is  just  as  good  as  new. 
This  situation  gives  rise  to  a  failure  probability  corresponding  to  the 
probability  of  disintegration  of  a  radioactive  nucleus.  The  element, 
like  the  nucleus,  has  no  memory  and  after  surviving  to  time  t  is  just  as 
good  as  it  was  at  time  0.  The  radioactive  decay  equation,  expressed 
in  symbols  familiar  to  the  reliability  specialist,  describes  the  time 
behavior  of  failure  of  the  element: 


-1  dN(t)  1  -t/t 

Cpi(t)  =  N(0)  dt  =  f «  (1) 

where  N(t)  is  the  number  of  units  that  have  not  failed  at  time  t  and  t 
is  the  mean  time  to  failure.  The  equation  is  in  the  normalized  form 
and  cp2(t)  is  expressed  in  fractional  probability.  The  functional  behavior 
of  <P2(t)  is  shown  in  Figure  1. 


If  the  element  is  used  a  comparatively  large  number  of  times  so 
that  it  begins  to  wear  out,  failure  of  a  different  type  occurs.  Just 
before  failure  in  this  mode,  the  element  is  not  as  good  as  new.  Its 
life  expectancy  is,  on  the  average,  quite  short  at  this  point  and  we 
may  expect  early  failure.  Failure  in  this  mode  will  scatter  about  a 
mean  in  such  a  way  as  to  approach  the  normal  or  Gaussian 
distribution,  which  expresses  the  distribution  of  random  fluctuations 
in  physical  phenomena.  The  functional  behavior  of  92^^^  shown  in 
Figure  2  and  is  given  by  the  equation: 


<P  2(t)  = 


tro 


(2) 


where  cr  is  the  dispersion  expressed  as  standard  deviation. 


If  the  chance  failure  curve  is  combined  with  the  wear-out  curve,  the 
function  of  Figure  3  is  obtained.  Here  the  shape  of  the  curve  is  similar 
to  the  density  function  expression  for  human  mortality,  which  is  also 
dotted  in  Figure  3.  The  shapes  of  the  two  curves  are  similar;  the 
major  difference  is  the  very  high  initial  value  of  the  human  mortality 
function.  This  initial  behavior  is  associated  with  accidents  at  birth, 
absence  of  immunity  to  disease  at  an  early  age,  and  with  abnormal 
weaknesses  of  the  very  yoiing,  and  has  no  counterpart  in  equipment 
being  discussed  here  unless  the  design  or  manvifacture  of  the  equipment 
is  fatilty. 

If  we  limit  our  discussion  to  short-lived  weapons  such  as  guided 
missiles,  wear  out  failures  are  in  general  negligible.  The  exceptions 
are  for  parts  subjected  to  severe  environmental  conditions.  An  example 
is  the  engine  of  the  V-1,  which  had  a  nominal  life  of  only  30  minutes 
because  of  severe  high  frequency  vibration.  Another  is  the  case  of  jet 
vanes  exposed  to  hot  exhaust  gases.  Still  another  exception  is  the  case 
of  the  absurd  "testing  to  death"  on  the  launching  pad  before  firing. 

We  are  then  left  with  chance  failure  as  the  primary  mode  of  failure 
of  significance  in  the  proper  operation  of  the  guided  weapon.  Is  it 
possible  to  elicit  additional  information  by  analysis? 

Stress  Levels 


If  we  examine  the  loading  of  an  element  such  as  a  beam  we  observe 
that  there  is  some  nominal  maximum  stress  to  which  the  beam  is 
ordinarily  subjected.  Here,  too,  the  element  of  chance  enters,  and  a 
distribution  about  whatever  maximum  stress  is  selected  will  surely  be 
observed.  The  normal  law  surely  does  not  exactly  fit  the  stress 
experience,  but  will  give  a  satisfactory  approximation  at  least  for 
semi -quantitative  discussion;  it  therefore  will  be  assumed  to  suitably 
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describe  the  scatter  about  the  nominal  maximxam.  That  the  stress  re¬ 
quired  to  cause  failure  of  the  element  will  similarly  vary  about  a 
nominal  value  is  well  known  to  those  who  have  concerned  themselves 
with  quality  control.  In  the  case  of  those  elements  for  which  a  failure 
strength  is  established  as  a  result  of  physical  testing,  the  maximum 
of  the  probability  function  will  be  located  at  a  higher  value  than  the 
nominal  because  strength  is  arbitrarily  defined  to  be  that  of  the  weakest 
member  of  the  group  of  elements  tested.  The  argument  here  will 
nevertheless  be  valid  with  appropriate  modification  because  this  pro¬ 
cedure  does  not  alter  the  existence  of  the  distribution  of  failure  strength, 
but  merely  adds  to  the  factor  of  safety. 


Figure  4  demonstrates  the  result  of  the  scatter  for  both  maximum 
stress  experience  expectancy  and  failure  strength  expectancy.  The 
nominal  maximum  stress  experienced  in  a  test  is  assumed  to  be  10 
units;  the  stressed  mexnber  is  designed  for  failure  at  a  nominal  20 
units,  giving  a  factor  of  safety  of  2.  We  observe  that  for  a  load  stress 
corresponding  to  f,  all  elements  characterized  by  a  failure  strength 
of  f  or  less  will  fail;  for  a  failure  strength  greater  than  f,  the  element 
will  survive.  Beyond  this  value,  the  stress  falls  short  of  causing  a 
failure.  Using  the  normal  density  function  (2),  the  failure  probability 
is  represented  by  the  expression: 
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where  a  is  the  stress  corresponding  to  the  nominal  maximum  stress 
ejqjerience;  b  the  stress  corresponding  to  the  nominal  failure  strength; 
and  a  ^  ando  2  the  standard  deviations  of  the  stress  experience  and 
failure  strength  curves,  respectively.  The  failure  probability  function 
Q(a,  b)  can  be  integrated  by  appropriate  change  of  variable,  with  the 
result  (3): 


00 


Q(a,  b)  = 


V  2  Tr(ai^+c^^)  b-a'^ 


]  df 


(4) 


This  integral  is  easily  evaluated  using  a  table  of  the  cumulative 
normal  distribution.  Figure  5  presents  the  graphical  interpretation 
of  Q(a,  b). 
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Stress-Strain  Variability  Relationships 


Examination  of  the  failure  probability  function  reveals  that,  once 
the  nominal  values  a  and  b  are  fixed,  it  is  a  function  of  the  standard 
deviations,  these  quantities  fixing  shape  and  therefore  the  extent  of 
Overlap  of  the  curves.  The  significance  of  this  statement  necessitates 
repetition.  Once  the  nominal  values  of  maxim\im  stress  experience  and 
of  failure  strength  are  fixed,  the  standard  deviations  or  ^  and  c  ^  deter¬ 
mine  the  value  of  the  failure  probability  function  Q(a,  b).  Furthermore, 
the  function  is  highly  sensitive  to  the  values  of  the  standard  deviations. 

Figures  4  and  6  illustrate  the  effect  of  the  magnitude  of  the  standard 
deviations  on  the  shape  of  the  curves.  Each  of  the  curves  is  normalized  ; 
that  is,  the  total  area  under  the  curves  is  unity  when  the  abcissa  is 
expressed  in  units  of  standard  deviation.  Li  the  example  represented 
by  the  curves  of  Figure  4  and  nominal  maximum  value  of  stress 
experience  is  a  =  10  stress  units;  the  nominal  value  of  failure  stress 
is  b  =  20  stress  units.  The  standard  deviations  cr  ,  and  (r  ^  ^-^e  each 
equal  to  4  stress  units.  The  value  of  the  function  Q(a,  b)  is  0.039, 
corresponding  to  a  reliability  of  96  per  cent.  In  the  example  represented 
by  the  curves  of  Figure  6,  nominal  stress  values  are  unchanged,  but 
standard  deviations  are  reduced  to  one  half  of  the  former  value,  e.  g. 

2  stress  units.  The  curves  are  much  more  sharply  peaked  about  the 
nominal  values.  Evaluation  of  the  failure  probability  function  for  this 
case  yields  a  failure  rate  of  2  per  10,  000,  down  by  a  factor  of  nearly  200 
from  that  of  the  first  example.  Further  reduction  of  each  of  the  standard 
deviations  to  one  stress  unit  vastly  improves  the  reliability,  the  failure 
rate  being  1  in  1  trillion.  If  the  density  functions  were  compressed  to 
the  limit,  i.  e. ,  a  delta  function  about  the  nominal  values,  the  failure 
probability  function  would  clearly  approach  zero  as  a  limit.  The  low 
value  obtained  for  Q(a,  b)  in  the  third  example  demonstrates  that  the 
failure  expectancy  becomes  vanishingly  small  long  before  delta  function 
liniit  is  reached. 

The  purpose  of  this  discussion  is  to  elucidate  the  effect  of  scatter 
on  the  failure  probability  function.  No  one  will  be  so  foolish  as  to 
purposely  design  a  part  with  a  failure  stress  less  than  the  expected 
stress.  The  failure  is  a  result  of  scatter.  Control  of  this  scatter,  or 
more  specifically,  the  magnitude  of  the  standard  deviations  is  the  key 
to  the  control  of  reliability  of  the  unit  under  consideration.  The  magni¬ 
tude  of  the  maximum  stress  experience  standard  deviation  may  be,  at 
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least  in  theory,  completely  at  our  disposal,  such  as  in  the  case  of 
screaming  in  a  rocket  engine;  it  may  be  completely  independent  of  our 
will,  as  in  the  case  of  countermeasures  against  our  weapon.  In  the 
first  case,  as  our  engineering  of  the  combustion  process  becomes  in¬ 
creasingly  precise,  we  reduce  the  standard  deviation  of  the  maximxim 
stress  experience  curve  by  being  able  to  more  closely  specify  just 
what  the  maximum  chamber  pressure  will  be.  In  the  second  case, 
about  all  that  can  be  done  is  to  specify,  for  example,  the  overpressure 
at  which  we  will  permit  an  airframe  structure  to  fail. 

The  magnitude  of  the  failure  strength  standard  deviation  is  com¬ 
pletely  at  our  disposal,  within  the  limits  of  our  design  and  production 
knowledge.  The  control  of  this  scatter  is  the  proper  function  of 
quality  control  activity.  It  is,  however,  a  matter  of  primary  impor¬ 
tance  that  the  quality  control  apply  to  the  actual  conditions  of  use  of 
the  item.  For  example,  the  scatter  of  failure  stress  for  a  non-scream¬ 
ing  rocket  engine  would  be  expected  to  be  less  than  that  of  one  that 
screams  because  of  the  greater  heat  transfer  coefficient  and  the  con¬ 
comitant  higher  wall  temperature  of  the  latter. 

Example 

To  summarize  the  discussion  of  the  scatter  of  the  maximvim  stress 
experience  and  failure  strength  functions,  we  may  state  that  good 
engineering  judgment  and  practice,  including  precise  quality  control, 
are  the  effective  means  to  control  of  reliability.  When  the  functions 
can  be  closely  specified,  the  factor  of  safety  can  be  reduced  to  a  low 
value  and  reliability  maintained.  Conversely,  when  scatter  of  these 
functions  is  large,  the  price  of  our  ignorance  is  an  increased  factor  of 
safety  or  a  reduced  reliability,  or  some  combination  of  the  two. 

At  this  point  a  mention  of  the  customarily  discussed  compound 
probability  relation  is  in  order.  The  over-all  reliability  of  an  assembly 
of  units  which  are  assvimed  to  be  mutually  independent  and  which  must 
all  function  to  avoid  failure  is  the  product  of  the  reliabilities  of  the 
individual  units.  Thus,  the  over-all  reliability  of  10  of  the  units  of 

99.  98  per  cent  individual  reliability  of  example  2  above  is  about  99.  8 
per  cent  and  still  better  than  the  single  unit  of  example  1.  Over-all 
reliability  of  100  units  is  98  per  cent,  and  of  1000,  82  per  cent.  For 
the  third  example,  the  compound  probability  of  failure  for  1000  units 
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is  only  about  1  in  1  billion.  The  importance  of  scatter  control  is  seen 
to  be  overwhelming;  by  reducing  the  standard  deviations  by  a  factor 
of  four,  failure  expectancy  is  reduced  by  a  factor  of  several  million. 
The  other  obvious  conclusion  is  that  reduction  of  the  number  of 
elements  of  an  assembly  is  of  exponential  weight  in  the  reliability 
problem.  Redvindancy  of  critical  components  is  an  extreme  measure 
that  can  be  used  to  reduce  the  failure  expectancy.  It  is  effectively 
applied  in  the  case  of  multi-engine  aircraft,  however. 


Summar'' 


To  summarize,  we  have  observed  that  the  failure  of  components 
of  guided  weapons  is  generally  in  the  nature  of  chance  rather  than 
wear-out  failure.  A  failure  probability  function  was  developed  to 
show  the  interrelation  of  maximum  stress  encountered,  component 
failure  strength,  and  the  scatter  of  the  values  of  these  parameters. 

It  was  shown  that  the  failure  probability  function  was  highly  sensitive 
to  scatter,  and  that  where  a  multiplicity  of  components  was  considered, 
a  moderate  reduction  of  scatter  could  result  in  the  difference  between 
success  and  failure.  The  same  result  can  be  achieved  by  reduction  of 
complexity. 
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DRi.'AL  DENSITY  FUNCTIONS  FOR  MAXIMUM 
IRESS  EXPERIENCE  AND  FAILURE  STRENGTH 


FIG.  5.-  FAILURE  PROBABILITY  FUNCTION  Q(o,b) 


EFFECT  OF  REDUCED  SCATTER  ON 
NORf.iAL  DENSITY  FUNCTIONS 


Appendix  A 

By  E.  G,  Pieruschka 

Classical  Theory  of  a  Combination 
of  Univariate  Stress  and  Strength  Distributions 

A  combination  of  univariate  stress-strength  distributions  is  shown  in 
figure  7 . 


^^st  ^  ^^th 

Figure?.  StreSs-Strength  Distribution. 


Here  we  have  a  stress  distribution  <}>gt  ( x)  with  the  mean  .  The  strength 
of  the  device  type  varies  with  device  to  device  accoi'ding  to  the  strength 
distribution  (x)  with  the  mean  p^h*  given  device  may  have  the  strength 
xj^.  The  stress  during  the  trial  of  the  device  may  be  Xg^.  When  xth  is 
larger  than  Xgt  we  have  a  successful  trial.  When  x^j^  is  smaller  than  Xg^  we 
have  a  malfunction.  Thus,  the  probability  of  malfunction  Q  may  be  obtained  by 
computing 

00  X 

Q  =  /  /  4>th(2)  <}>st(5c)  dzdx 

x=—  00  X  —  -00 


(1.1) 


The  integration  area  of  the  two  dimension  function  <t>g^  (x)  (z)  is  indicated 


We  recognize  from  figure  8  that  the  double  integral  (1.1)  which  has  one  variable 
integral  limit  may  be  transformed  in  a  double  integral  with  fixed  integral 
limits: 
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Q=  /  /  (‘•2> 

v=-oo  u=  -00  \ 

When  <|)g^  (x)  and  (z)  are  normal  distributions  we  may  write  them  as: 


exp 


^  o-gt  I  2crgt^  J 

<t>th(2)  =  —  exp  [-  (z  -  pth)^  1 

o’th  L  2(rth  J 

Putting  (1.3)  and  (1.4)  into  (1.2)  we  have 


(1.3) 

(1.4) 


Q  = 
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^nij-gtO-th 


/  /  expf-n(HpL_. 

=  00  u=  00  I  ^  ^’st  O'st  '  'V^o-th  <nh  '  J  J 

(1.5) 


The  qiiantities  in  parentheses  may  be  written  as 

/  u-v  -  p-aL\  ,  respectively. 

\  ^  o^st  /  '■  ' 

In  order  to  simplify  the  left-hand  parenthesis,  we  substitute 

u=Vi  °'st  y  ^  M-st 

du  =V2  o-gt  dy 


Equation  (1.5)  becomes 
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On  simplifying,  the  exponent  becomes 

.  i  li-fyZ  +  ZylSl  ^(2v-V2d)^] 
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Equation  (1.9)  becomes 
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Since  for  any  value  of  v, 
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Equation  (1.  10)  becomes 
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where  or  =  ^(T  th^  +  o’  st'^ 


(1.10) 


(1.11) 


(1.  12) 


Due  to  symmetry  of  the  normal  distribution,  the  integration  limits  may  be 
reversed  with  opposite  sign,  giving  the  same  result  as  equation  (1.  3). 


Appendix  B* 


By  H.  A,  David 
Professor  of  Statistics 
Virginia  Polytechnic  Institute 


Probabilistic  Derivation  of  Probability  that  Stress  Exceeds  Strength 

Suppose  that  x,  denoting  stress,  is  normally  distributed  with  mean 
a  and  standard  deviation  and  that  y,  denoting  strength,  is  normally 

distributed  with  mean  b  and  standard  deviation  0“2*  The  probability 
of  failure  when  stress  exceeds  strength  is: 


Q  =  Pr  (x  >  y) 

=Pr{x-y-a+b> 


=  Pr  [.x_-.X-,-__a._-Lb.  > 
+  02^) 


=  Pr  (u  > 


b  -  a  • 


+  02^) 

where  u  is  unit  normal. 


-a  +  b) 

-  a  +  b  2 

■  )  , 


(2.1) 
(2.  2) 

(2.3) 
(2.  4) 


♦  This  appendix  gives  the  same  result  as  that  of  Appendix  A,  but  through 
a  different  approach. 
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